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Abstract 


Morbidity is a diseased state or symptom in which a person is unable to do daily activities. 
The morbidity rates can differ by person and place. This paper aims to investigate factors 
affecting morbidity rates in each province in Indonesia. Using applied spatial error model 
(SEM), this study analyzed the spatial effect on factors being investigated. It demonstrated 
that SEM is well-applied on Indonesian morbidity rates data. The results show that in 2020, 
provinces on Java Island and some provinces on Sulawesi Island have health problems that 
need to be addressed because they are in the high and moderately high morbidity categories. 
The value of the Moran Index is positive, indicating a similarity in the percentage of 
morbidity rates in adjacent provinces. The morbidity rates in most of provinces in Indonesia 
are affected by the duration averages of school attendance, health insurance, and population 
density. The coefficient of determination (R?) of the SEM model of 77.60% can be said to be 
a fairly good model. 


Keywords: morbidity, spatial error model, school attendance, health insurance, provincial minimum wage 
system 


Article History: 
Received: May 8, 2023 Revised: June 9, 2023 
Accepted: June 17, 2023 Published online: July 31, 2023 
Suggested Citation: 


Irma Yahya, Bahriddin Abapihi, Makkulau, Agusrawati & Gusti Ngurah Adhi Wibawa (2023). Spatial Error 
Model to Analyze Morbidity Rate in Indonesia. International Journal of Science, Technology, Engineering and 
Mathematics, 3 (3), 18-35. https://doi.org/10.53378/352998 


About the authors: 
‘Bachelor of Science in Statistics, Lecturer, Halu Oleo University, Indonesia. 


Corresponding author. Bachelor of Science in Statistics, Lecturer, Halu Oleo University, Indonesia. 


Corresponding email: gusti.a.wibawa(@uho.ac.id 


© The author (s). Published by Institute of Industry and Academic Research Incorporated. 
©) @ This is an open-access article published under the Creative Commons Attribution (CC BY 4.0) 
MEM Cense, which grants anyone to reproduce, redistribute and transform, commercially or non- 
commercially, with proper attribution. Read full license details here: 
https://creativecommons.org/licenses/by/4.0/. 


ISSN 2799-1601 (Print) 2799-161X (Online) | 19 


1. Introduction 

The level of public health can be described by morbidity, mortality, and nutritional 
status (Ministry of Health RI, 2010). According to Tulchinsky and Varavikova (2014), 
morbidity (measure of disease occurrence) can be measured based on the number of people 
with disease, period and duration of illness, frequency of deaths, disease, disability, and 
risk factors related to health outcomes. On the other hand, mortality or death rate is the 
measure of death frequency over a given period (Nolte & McKee, 2012) while nutritional 
status refer to the physiological state of the population in reference to nutrient intake (Food 
Agric Organ, 2007). One of many indicators to assess successful development in public 
health is the level of morbidity rates in the communities. The lower the morbidity, the 


healthier the people are. 


Morbidity rate is the number of certain disease incidences that is formulated as the 
number of people suffering from disease per 1000 people exposed to the disease (Kardjati 
et al., 1985; Singh-Manoux et al., 2008). The morbidity rate is more important to tackle 
than mortality rate, because high morbidity rate can also trigger high mortality rate. 
Morbidity rate can be used to measure the health condition in general, to assess how 
successful the programs of disease eradication, sanitation of environment, and to 


understand the people knowledge on health services. 


Based on research conducted by health experts (Okoroiwu, 2020; Moise, 2018; 
Wunsch & Gourbin, 2018; Zylke & Bauchner, 2020; Amini et al., 2021; Chang et al., 
2022), morbidity is caused by neonatal respiratory distress syndrome, tuberculosis and 
diarrhea. In most cases, asthma, tuberculosis and diarrhea diseases have a negative impact 
on the patient's life, causing children to often miss school, limiting personal and family 
activities and decreasing work productivity. The higher the morbidity rate, the worse the 
level of public health. The morbidity rate can reflect the real state of health because it has a 
close relationship with environmental factors such as poverty, malnutrition, infectious 
diseases, housing, proper drinking water, environmental hygiene and health services 
(Kardjati, et.all. 1985; Krieger & Higgins, 2002; Boyles et al., 2021; Jian et al., 2017; 
Kumari et al., 2023). 
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Based on data from Bangdan Pusat Statistik (BPS), the percentage of Indonesia's 
morbidity rate from 2017 was 14.31%, decreased to 13.19% in 2018 but experienced an 
increase again in 2019 to 15.38%. According to Soleman (2020), every five children in 
preterm birth in neonatal phase carried co-morbid factors that increase morbidity rate. The 
study of Hussain et al. (2015) found that one-third of the Indonesian adult population have 
multimorbidity. Aside from women being particularly affected, the study showed high 
prevalence of multimorbidity among younger individuals. Based on these data, the 
government must take steps to deal with it so that the morbidity rate does not increase in 


the following years. 


Morbidity rate differ by geographical area and country. It depends on life quality of 
the people within the area. The factors affecting morbidity rate also differ by geographical 
area. Accordingly, the method to analyze such a situation should incorporate spatial term 
in its model. One of the models that can tackle this condition is a spatial regression model. 
The spatial effect can show the clustering effects on adjacent areas. By this method, the 
spatial effects can be seen and be applied to analyze the morbidity rate in Indonesia using 


dataset obtained in 2020. 
2. Literature review 


2.1. Spatial Regression Models 

According to Rey et al. (2020), regression and prediction examines how spatial 
structure analyze data. Through spatial structure, regression models generate explicitly 
spatial data. If the model systematically mispredicts, a better model can be developed. For 
example, mapping classification or prediction error can help errors in data clusters. Hence, 
“regardless of whether or not the true process is explicitly geographic, additional 
information about the spatial relationships between observations can make predictions 


better” (Rey et al., 2020). 
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Spatial Dependence Test. According to Anselin (1988), the tests to be used to 
investigate the spatial dependence in the error term of a model are Moran’s J and Lagrange 
Multiplier (LM). Moran’s index is a kind of correlation to investigate the relationship 
among adjacent observations. Under null hypothesis the observations are assumed to have 
no spatial correlation (Ho: J = 0). The test statistic used is as follows. 

I-E(1) 
Z itun, So 
"me var) (1) 
n(n? —3n+3)S, —nS, +28.) 


i 
eas a R RE EE 


with = 
Accordingly, Moran’s index (Moran’s /) can be calculated as: 
Sod)” 
(2) 


n n * 1 n x s2 k x 7 
With a W Boa A Wy EW S=), (W +W 


i=1 Loy j=l j 


Where: 
y; : i-th observation (i= 1,2, ..., n) 
y; : j-th observation (j= 1, 2, ... n) 
y : general mean of observations 
W; : spatial weight matrix 
n : number of areas 
I : value of Moran’s I 
EV) : mean value of J 
var(/) : variance of I 


Reject Ho if the value of calculated Z greater than tabulated Zw/2. 


Spatial Heterogeneity Test. To test the effect of spatial heterogeneity, it can be done 
by applying Breusch-Pagantest (BP test). The hypotheses are: 
Ho: No heterogeneity among areas 


Hı : Heterogeneity exists among areas 
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The formula for Breusch-Pagan test is: 
se 1 T Tryn- gT 
BP = a ZZ Z) Zf 8) 


Where e7 is errorof i-th observation and Z is a matrix of independent variables of size n x (k 
p 


+1). Reject Ho if the value of BP>%?%a p. 


Lagrange Multiplier (LM) Test. LM test is applied to determine the existence of 
spatial effect in the model. The procedures of LM test are as follows. 
a) Lagrange Multiplier Lag 

The hypotheses are: 

Ho: p =0 (no spatial-lag dependency) 

Hi: p #0 (spatial-lag dependency exists) 

Test-statistic: 

Mo aa EN 
s (CWXB)’ MCWXB) + Ts~) 

(4) 
with 
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Decision criterion, reject HO if LM pg > i or p-value < a. 


b) Lagrange Multiplier Error 
The hypotheses are: 
Ho: å =0 (no spatial-error dependency) 
Hı : å #0 (spatial-error dependency exists) 
Test-statistic: 
(E Wels’? 
a T (5) 
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úl . . . . ph 
Decision criterion, reject HO if LM por > Xia X p-value < a. 


2.2. Spatial Regression 

Spatial regression is a regression in which it incorporates spatial effect which is 
represented in the spatial weight matrix based on areal adjacency. 

Spatial Model. Model for spatial regression according to Anselin is as follows: 


y=pWy+Xß+u 


(6) 
u=4Wu+e 
(7) 

Substituting equation (6) to (7) becomes 
y=pWy+XBp+AWute e€~N(0,071 (8) 
Where, 

Y __: vector for response variable (nx/) 

X  : matrix for predictor variables (nx(k+J)) 

B —_: coefficient vector of regression parameters 

P  : Spatial lag coefficient parameter of dependent variable 

A  : Spatial lag coefficient parameter on errors 

u : vector for error terms on Y (n x 1) 

n  : number of observations or locations 

I : identity matrix (n x n) 

g : vector for error term in equation U which has normal distribution with 


mean zero and variance o’I of size n x 1 
Spatial regression with area approximation consists of some models i.e., Spatial Error 


Model(SEM). General model for SEM can be seen as follows: 


y=XB+4Wu+s e~N(0,071) (9) 


Where AWu shows spatial on spatially dependent error(s). 
Parameter estimation is obtained by maximizing the logarithm or using maximum likelihood 


method. Parameter estimation of SEM as the following equation: 
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B = (X"(I-2W)' (I- AW)X)|(X"(I- AW)’ (I- 2W)y) (10) 


2.3. Parameter Significance of Spatial Regression 


Model Fitness Test. To test the fitness of SAR model, the following procedure is 
employed: 
Ho: p= 8, =0 
Hı : o #0 (at least there isa 2, #0 for k= 1, 2,..., p) 


Test statistic used is: 


$ JKRIk | R° Ik 
mwe JKĶGI/(n-k-1) (1-R°)/(n-k-1) 


(1) 

Where: 
k : number of predictor variables 
R? : coefficient of determination 


Decision criterion: reject Ho if Featculated> F(a,k,n-k-1). 


Partial Test. The Hypothesis for partial test is as follows: (Anselin,2003) 
H,:@9=0 (Parameter is unsignificant) 
H,:@0#0 (Parameter is significant) 
Test statistic for parameter significance is: 
6 


Z nitung = sbo (12) 


Where Sb, D be asymptotic standard error, and @ be the parameter of spatial regression 


(namely $, 2 and p). If Znitung> Zoor p-value<a (a =0,05), then the decision is to reject 


Ho, which means the coefficient of regression can be used in the model. 
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2.4. Selecting the Best Model 
The criterion of model selection, which is used in this paper, is AIC criterion. The 
model is selected when it has the least AIC score. As presented by Hu (2007), the formula to 


calculate the AIC score is as follows: 


AIC = -2log(L(6| y)) +28 
(13) 


Where: 
L(6| y) : likelihood function of the parameter being estimated. 


2.5. Morbidity Rate 


Morbidity is the condition of a person who is said to be sick if the health complaints 
that are felt cause disruption of daily activities, namely not being able to work, take care of 
the household, and normal activities as usual. According to Hernandez and Kim (2022), 
morbidity is one of the two commonly measures of epidemiological surveillance which 
describes the progression of a given health occurrence. Morbidity is usually represented or 
estimated using prevalence (proportion of population) or incidence (frequency within a 


population). It can be presented as a ratio or percentage (Mont, 2014). 


The following is the morbidity rate formula. 


aMm=-KK x100 
JP 
(14) 
Description: 


AM : morbidity rate 
JPKK : number of people who experience health complaints and disruption of activities 
JP : total population 

According to Choi et al. (2019), morbidity indicators include prevalence, incidence 
and attack rate. Prevalence measures the proportion of population with specific illness at a 
specific time period which can be point or period prevalence. Point prevalence measures the 
number of cases in the population at a point in time while period prevalence estimates the 


proportion of individuals with a certain condition at any time during a specified time period. 
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Hernandez and Kim (2022) refers incidence to the frequency at which individuals within a 
specific population develop a given symptom or quality. It could account the number of 
people with newly developed disease (Choi et al., 2019). Meanwhile, attack rate is the 
proportion of incident cases of a disease occurring to population exposed to the source of 
disease. 

There are three dimensions that show indicators of morbidity, namely the long and 
healthy life dimension, the knowledge dimension and the decent life dimension. The 
longevity and health dimension is measured based on life expectancy and the percentage of 
population seeking medical treatment at a health worker's practice. The knowledge 
dimension is measured based on the illiteracy rate and the average years of schooling of the 
population 15 years and above. While the dimension of a decent life is measured based on 
the percentage of the population that has access to decent drinking water and the percentage 
of the poor population. 

The Pan American Health Organization (n.d.) enumerates the factors affecting the 
accuracy of the measurement of morbid events such as data quality, validity of measurement 
instruments, disease severity, cultural norms, confidentiality and health information system. 
It is particularly important that the analysis take into consideration the diversity and volume 
of data. Since data differ by area of a country, it is imperative that the model used identifies 
the disaggregation of data. Similarly, the probability of errors in the data analysis must be 
treated with utmost confidence. 

2.6. Social Factors 

Social factors are factors that are born, grow, and develop in a common life (Salim, 
2002). According to Anderson (1995), social factors include education and ethnicity. 
Increasing education is an effort to prevent morbidity. In a research conducted by 
Ardhiyanti (2013), the level of education in an area is reflected in the illiteracy rate 
because it can show the inability of the population in an area to absorb information from 
various media, communicate orally and in writing. In addition to the illiteracy rate, in 2014 


BPS used the average years of schooling as an indicator of education. 


According to Institute of Medicine (2001), socioeconomic status is a strong and 
consistent predictor of morbidity. In fact, most studies showed that lower morbidity rate 


occurs among socioeconomically advantaged people (Kaplan & Keil, 1993; Syme & 
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Berkman, 1976; Lantz et al., 1998; Caldwell, 1990; Lantz el al., 2010; Bosworth, 2018; 
Dhrifi, 2018). 


2.7. Environmental and Behavioral Factors 


Morbidity reflects the actual state of health because it has a close relationship with 
environmental factors such as malnutrition, infectious diseases, housing, healthy drinking 
water, environmental hygiene, and health services (Wulandari, 2017). According to Pan 
American Health Organization (n.d.), there are environmental factors that are difficult to 
measure such as exposure to air pollution, exposure to sunlight and exposure to some disease. 


3. Methodology 
Data used in this research were obtained from Allstats application of BPS (2020). The 
variables were morbidity rates (Y), average duration of school attendance (X1), proper 
drinking water availability (X2), proper sanitation (X3), health insurance (X4), population 
density (Xs), minimum provincial wage (X6). 
The research was conducted based on the following steps. 
e Summarizing descriptive statistics of variables 
e Analyzing SEM model as follows: 
o Performing spatial effect test and model identification 
o Conducting assumption test of regression for the following SEM model 
o Interpreting model and making conclusion 
There are several studies that model mordibitas cases, as follows: 
Factors Affecting Morbidity of East Java Residents with Multivariate Geographically 
Weighted Regression (MGWR) (Hanum, 2013). 
Modeling Factors Affecting Morbidity in Central Java Using Spline Truncated 
Nonparametric Regression (Rosanti, 2020) 


4. Findings and Discussion 


4.1. Description of Morbidity Rate 
For simplicity in interpreting the spread of morbidity in each province in Indonesia, 
we do the mapping as seen in Figure 1, where brighter color represents high morbidity rate, 


while lighter color shows the opposite. 
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Figure 1 
Map of Morbidity Spread in Indonesia 2020 


Y 


[7.810 : 11.340] (8) 
E [11.520 : 13.550] (9) 
E [13.610 : 14.710] (9) 
E [15.050 : 21.090] (8) 


Based on Figure 1, it can be seen that there are 8 provinces to have very high 
morbidity rates and all these provinces are located in Java. Some 8 other provinces have 
high morbidity rates. This means that provinces on Java Island and some provinces on 
Sulawesi Island have health problems that need to be addressed because they are in the 
high and moderately high morbidity categories. There are 9 provinces of moderate 
morbidity rates which are distributed over Sumatera, Kalimantan, and Sulawesi islands, 


and the other 8 provinces with low morbidity rates. 


4.2. Spatial Effect Tests 
In SEM model, there are two spatial effect tests to be stratified, namely spatial 
dependence and spatial heterogeneity tests. The results of the tests are presented on Tables 


1 and 2. 


Table 1 
The Result of Moran Index Test 


Morans’I p-value Pattern 


0.357063 0.017459 Grouping 
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Based on table 1, the p-value is 0.017459 < a which means H, rejection and there 


is spatial autocorrelation among adjacent provinces. The Moran’s index is 0.357063 which 
shows positive autocorrelation. This means there are similarities of morbidity rates among 


adjacent locations and tends to grouping. 


Table 2 
The Result of Spatial Heterogeneity Test 


Breusch-Pagan Probability Description 


3.7045 0.71659 Accep H, 


Based on table 2, the probability of the BP test for this model is 0.71659 which is 
greater than 0.05, so no rejection of Ho. This result means that no heterogeneity among 


locations. 


4.3. Initial Identification of Spatial Model 
LM test detects the existence of spatial dependence more specific in term of lag, 
error, or lag and error. If LM (lag) is significant then the model is Spatial Autoregressive 
(SAR). But, if LM (error) is significant then the model is Spatial Error (SEM). While if 
both are significant then the model is Spatial Autoregressive Moving Average (SARMA). 


Table 3 
Result of LM Tests 
Uji LM Statistic p-value Description 
Lagrange Multiplier (lag) 1.1780 0.27776 Accept H, 
Lagrange Multiplier (error) 4.9320 0.02636 Reject H, 


Based on table 3, p-value of LM (error) test is less than 0.05 then the decision is to 
reject H,. This means there is spatial error dependence on morbidity and we continue to 


SEM model construction. 
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4.4, Spatial Error Model (SEM) 
In model detection, it was concluded that the appropriate model is SEM. Hence, we 
are going to use SEM. The following is the result of SEM output based on parameter 


estimation with level of significance at 5%. 


Table 4 
Parameter estimation dan testing of SEM model 
Parameter Coefficient Probability 
Constant 32.8322 0.0000 
Xi -1.8299 0.0000 
X2 - 0.0754 0.0018 
X5 0.0004 0.0007 
Lambda 0.6601 0.0000 


Based on table 4, it can be seen that the variables affecting morbidity in Indonesia 
are average duration of school attendance (X1), health insurance (X4) and population 


density (Xs). Using these significant variables, SEM model is then constructed. 
The equation of SEM model is as follow. 


Y = 32,8322 — 1,8299X, — 0,0754X, + 0,0004X; + U; 
34 
dengan: U; = 0,6601 ` W;;jUj; + £i 
j=Lj+i 
Based on the model with some significant variables, the interpretations are as 


follows: 


1. The coefficient for average duration of school attendance is -1.8299 which means if 
we increase the average duration of school attendance for one year and other 
variables remain constant, then the morbidity rate in Indonesia will decrease 1.8%. 

2. The coefficient for health insurance is -0.0754 which means if this variable increase 
ten percent, the morbidity rate tends to decrease 0.75%. 

3. The coefficient for population density is 0.0004 which means if the population 


increases 100% then the morbidity rate will also increase as much as 0.04%. 
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Coefficient of determination (R?) of the SEM model is 77.60 % which represents 


that the model is good enough to explain the morbidity rate in Indonesia. 


5. Conclusion 


Based on the Moran’s index testing, it can be concluded that there is spatial 
autocorrelation on morbidity rates, which means similarities among adjacent provinces 
exist. The variables that affect morbidity rates in Indonesia are average time duration of 
school attendance (X1), health insurance (X4) and population density (Xs), where the SEM 


model can be written as: 


A 


Y = 32,8322 — 1,8299X, — 0,0754X, + 0,0004X. + Ui 
34 
with U; = 0,6601 > W;;jU; + £i 
j=1Lj+i 
From the model obtained, it can be seen that the factors significantly affect morbidity 
cases in Indonesia. The results imply the need for the government to further formulate 
measures in reducing the percentage of morbidity, especially the factors that have the greatest 
influence, namely the average length of schooling and health insurance. The results clearly 
show that strengthening the education of the population greatly reduces morbidity rate. 
Taking this in the context, it is a challenge to the government to further improve the quality 
of education and health insurance in the country. The model suggests a country with healthy 
mind and healthy body to reduce morbidity rate. Given the limits of the study, further 
research is recommended considering the number of predictor variables and types of 


variables to obtain a better model. 


6. Acknowledgement 
Authors would like to thank the Ministry of Research and Technology and Halu 


Oleo University for funding this research project from 2021. 


References 


Amini, M., Zayeri, F. & Salehi, M. (2021). Trend analysis of cardiovascular disease 


mortality, incidence, and mortality-to-incidence ratio: results from global burden of 


32 | International Journal of Science, Technology, Engineering and Mathematics, Volume 3 Issue 3 


disease study 2017. BMC Public Health 21, 401 (2021). 
https://doi.org/10.1186/s12889-021-10429-0 

Anselin, L. 1988. Spatial Econometric: Methods And Models. Dordrecht: Kluwer 
Academic Publishers. 

Bosworth, B. (2018). Increasing Disparities in Mortality by Socioeconomic Status. Annual 
Review of Public Health 2018 39:1, 237-251 

Boyles, A. L., Beverly, B. E., Fenton, S. E., Jackson, C. L., Jukic, A. M. Z., Sutherland, V. 
L., Baird, D. D., Collman, G. W., Dixon, D., Ferguson, K. K., Hall, J. E., Martin, E. 
M., Schug, T. T., White, A. J., & Chandler, K. J. (2021). Environmental Factors 
Involved in Maternal Morbidity and Mortality. Journal of women's health (2002), 
30(2), 245-252. https://doi.org/10.1089/jwh.2020.8855 

BPS (2020). Angka Harapan Hidup (AHH) MenurutProvinsi dan JenisKelamin. 
AplikasiAllstats, Indonsia. 

BPS (2020). Angka Morbiditas Indonesia 2020. Aplikasi Allstats, Indonesia. 

BPS (2020). Persentase Penduduk Miskin Menurut Provinsi. Aplikasi Allstats, Indonesia. 

BPS (2020). Rata-Rata Lama Sekolah. Aplikasi Allstats, Indonesia. 

BPS (2020). Rumah Tangga Menurut Provinsi, Tipe Daerah dan Sumber Air Minum Layak. 
AplikasiAllstats, Indonesia. 


Caldwell, J. C. (1990). Cultural and Social Factors Influencing Mortality Levels in 
Developing Countries. The ANNALS of the American Academy of Political and Social 
Science, 510(1), 44-59. https://doi.org/10.1177/00027162905 10001004 

Chang, D., Chang, X., & He, Y. (2022). The determinants of COVID-19 morbidity and 
mortality across countries. Sci Rep 12, 5888 (2022). https://doi.org/10.1038/s41598- 
022-09783-9 

Choi, J., Ki, M., Kwon, H. J., Park, B., Bae, S., Oh, C. M., Chun, B. C., Oh, G. J., Lee, Y. H., 
Lee, T. Y., Cheong, H. K., Choi, B. Y., Park, J. H., & Park, S. K. (2019). Health 


Indicators Related to Disease, Death, and Reproduction. Journal of preventive 
medicine and public health = Yebang Uihakhoe chi,52(1), 14-20. 
https://doi.org/10.396 1/jpmph. 18.250 

Committee on Population; Division of Behavioral and Social Sciences and Education; Board 
on Population Health and Public Health Practice; Health and Medicine Division; 


National Academies of Sciences, Engineering, and Medicine. Improving the Health 


ISSN 2799-1601 (Print) 2799-161X (Online) | 33 


of Women in the United States: Workshop Summary. Washington (DC): National 
Academies Press (US); 2016 Mar 23. 3, Socioeconomic and Behavioral Factors That 
Influence Differences in Morbidity and Mortality. Available from: 
https://www.ncbi.nlm.nih.gov/books/NBK356223/ 


Dhrif, Abdelhafidh (2018). Health-care expenditures, economic growth and infant mortality: 
evidence from developed and developing countries. CEPAL Review No. 125 

Food Agric Organ. Nutritional Status Food Security. Food Agric Organ. (2007). Available 
online at: https://elearning.fao.org/course/view.php?id=189 


Hanum, D., & Puhardi (2013). Faktor-Faktor Yang Mempengaruhi Morbiditas Penduduk 


Jawa Timur Dengan Multivariate Geographically Weighted Regression (MGWR). 
Jurnal Sains dan Seni POMITS. Vol. 2, No. 2. 

Hernandez, JBR & Kim, PY (2022). Epidemiology Morbidity And Mortality. In: StatPearls 
[Internet]. Treasure Island (FL): StatPearls Publishing. Available from: 
https://www.ncbi.nlm.nih.gov/books/NBK547668/ 

Hussain MA, Huxley RR, & Al Mamun A. (2015). Multimorbidity prevalence and pattern in 
Indonesian adults: an exploratory study using national survey data. BMJ Open 
2015;5:e009810. doi: 10.1136/bmjopen-2015-009810 

Institute of Medicine (US) Committee on Health and Behavior: Research, Practice, and 
Policy. Health and Behavior: The Interplay of Biological, Behavioral, and Societal 
Influences. Washington (DC): National Academies Press (US); 2001. 4, Social Risk 
Factors. Available from: https://www.ncbi.nlm.nih.gov/books/NBK43750/ 

Jian, Y., Messer, L. C., Jagai, J. S., Rappazzo, K. M., Gray, C. L., Grabich, S. C., & Lobdell, 
D. T. (2017). Associations between Environmental Quality and Mortality in the 
Contiguous United States, 2000-2005. Environmental health perspectives, 125(3), 
355-362. https://doi.org/10.1289/EHP119 

Kardjati S., Alisjahbana A, A. & Kursin J. (1985). Aspek Kesehatan dan Gizi Anak Balita. 
Akarta: Yayasan Obor Indonesia. 

Kemenkes, RI (2010). Profil Kesehatan Indonesia 2009. Kemnterian Kesehatan Republik, 


Jakarta 2010. 
Krieger, J., & Higgins, D. L. (2002). Housing and health: time again for public health action. 
American journal of public health, 92(5), 758-768. 


https://doi.org/10.2105/ajph.92.5.758 


34 | International Journal of Science, Technology, Engineering and Mathematics, Volume 3 Issue 3 


Kumari U, Sharma R, & Keshari J. (2023) Environmental Exposure: Effect on Maternal 
Morbidity and Mortality and Neonatal Health. Cureus 15(5): e38548. 
doi: 10.7759/cureus.38548 

Lantz PM, House JS, Lepkowski JM, Williams DR,Mero RP,Chen J. (1998). 
Socioeconomic Factors, Health Behaviors, and Mortality: Results From a Nationally 
Representative Prospective Study of US Adults. JAMA. 1998;279(21):1703—1708. 
doi: 10.1001/jama.279.21.1703 

Lantz, P. M., Golberstein, E., House, J. S., & Morenoff, J. (2010). Socioeconomic and 
behavioral risk factors for mortality in a national 19-year prospective study of U.S. 
adults. Social science & medicine (1982), 70(10), 1558-1566. 
https://doi.org/10.1016/j.socscimed.2010.02.003 

Lesage, J.P. (1999). The Theory And Practice Of Spasial Econometrics. Departement of 


Economics University Of Toledo. 

Moise, I. (2018). Causes of Morbidity and Mortality among Neonates and Children in Post- 
Conflict Burundi: A Cross-Sectional Retrospective Study. Children, 5(9), 125. MDPI 
AG. Retrieved from http://dx.doi.org/10.3390/children5090125 

Mont, D. (2014). Morbidity Measures. In: Michalos, A.C. (eds) Encyclopedia of Quality of 


Life and Well-Being Research. Springer, Dordrecht. https://doi.org/10.1007/978-94- 
007-0753-5_ 1849 
Nolte E. & McKee M. (2012). In amenable mortality — deaths avoidable through health care 


— progress in the US lags that of three European countries. Health Aff. 2012;31:2114— 
2122. 

Okoroiwu, H. U., Uchendu, K. I., & Essien, R. A. (2020). Causes of morbidity and mortality 
among patients admitted in a tertiary hospital in southern Nigeria: A 6-year 
evaluation. PloS one, 15(8), e0237313. https://doi.org/10.1371/journal.pone.0237313 

Pan American Health Organization (n.d.). Health Indicators: Conceptual and Operational 
Considerations. A Publication of the World Health Organization Regional Office for 
the Americas 


Rey, S.J., Arribas-Bel, D. & Wolf, L.J. (2020). Geographic Data Science with Python. 
Available online at https://geographicdata.science/book/intro.html# 


ISSN 2799-1601 (Print) 2799-161X (Online) | 35 


Rosanti, I, W., & Budiantara, I, N. (2020). Pemodelan Faktor-Faktor Yang Mempengaruhi 
Morbiditas di Jawa Tengah Menggunakan Regresi Nonparametrik Spline Truncated. 
Jurnal Inferensi. Vol. 3, No. 2. 

Sani Rachman Soleman (2020). The Trend of Children Mortality Rate In Indonesia. Jurnal 
Ilmu Kesehatan Masyarakat, Mar 2020, 11(1):52-62. 
https://doi.org/10.26553/jikm.2020.11.1.52-62 

Singh-Manoux A, Gue'guen A, Ferrie J, Shipley M, & Martikainen P. (2008) Gender 
Differences in the Association Between Morbidity and Mortality Among Middle- 
Aged Men and Women. Am J Public Health 98: 2251—2257. 

Tulchinsky, T. H., & Varavikova, E. A. (2014). Measuring, Monitoring, and Evaluating the 
Health of a Population. The New Public Health, 91-147. 
https://doi.org/10.1016/B978-0-12-415766-8.00003-3 

Wulandari, K. 2017. Pemodelan Faktor-Faktor Yang Mempengaruhi Angka Morbiditas di 


Jawa Timur Menggunakan Regresi Nonparametrik Spline [Tesis]. Surabaya: Institut 
Teknologi Sepuluh Nopember, Fakultas Matematika dan Ilmu Pengetahuan Alam. 
Wunsch, G. & Gourbin, C. (2018). Mortality, morbidity and health in developed societies: a 
review of data sources. Genus 74, 2 (2018). https://doi.org/10.1186/s41118-018- 
0027-9 
Zylke JW. & Bauchner H. (2020). Mortality and Morbidity: The Measure of a Pandemic. 
JAMA. 2020;324(5):458—459. doi:10.1001/jama.2020.11761 


